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Abstract 



This paper presents a set of quantum Reed-Muller codes which are typi- 
cally 100 times more effective than existing quantum Reed-Muller codes. The 
code parameters are [[n,k,d]} = [[2 m , £[ =0 C(m, I) - E^/" 1 C(m, Z), 2 m " r ]] 
where 2r + 1 > m > r. 
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Quantum information processing, which includes quantum communication, cryp- 
tography, and computation is currently moving from theoretical analysis towards 
implementation [|TJ . The fragility of quantum states has led to the the development 
of quantum error correcting codes || @ • A large number of research papers have ap- 



ed 

peared since it was shown that quantum error correcting codes exist P|-|l2|. Much 
of this research work has focused on repetition codes, i.e., codes with parameters 
[[n,l,d]], or single error correcting codes, i.e., codes with parameters [[n,k,3]]. Rep- 
etition codes are inefficient, with code rates (R = -). Single error correcting codes 
are not very effective because they correct only one error in a block of n qubits. It 
is desirable to design quantum error correcting codes that are more efficient, and 
more powerful. 

Calderbank and Shor described a general method to construct non-trivial, multiple 
error correcting quantum codes || . The same method was independently discovered 
by Steane |J. Using this method, together with a technique to accommodate more 
information qubit(s) in a coded block by slightly reducing the minimum Hamming 
distance, Steane presented a set of multiple error correcting codes based on trivial 
codes and Hamming codes [|10j, and a set of quantum Reed-Muller Codes based on 
classical Reed-Muller codes [BJ. In this paper, a new set of quantum Reed-Muller 
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codes are derived using the technique of Calderbank & Shor/Steane. These codes 
have lower code rates than Steane's but their minimum Hamming distances are 
larger. It will be shown that if the uncoded quantum system has a qubit error rate 
less than 0.3%, one particular code of rate R = 0.246 will be able to bring the output 
qubit error rate down to 10~ 9 . 



2. Quantum Codes 

A multiple error correcting quantum code described in Q comprises two classical 
error correcting codes, C\ = (n,ki,di) and C 2 = (n,k 2 ,d 2 ). C\ and C 2 are related 
by 

{0} C C 2 C C x C F 2 n (1) 

where {0} is the all-zero code word of length n, and F 2 is the n- dimensional binary 
vector space. The quantum code C defined by G\ and C 2 has the parameters 

[[n, k, d\] = [[n, dim(C^) - dim(C 1 ± ), min(c/ 1 , d 2 )}} (2) 

where C 1 - represents the dual code of C; the parameters of C 1 - are (n,n — k^d 1 -); 



d and d are indirectly related by Mac Williams theorem p3| ; and dim(C) denotes 
the dimension of the code C. 

In order for the error correction scheme to work, coded quantum states must be 
represented in two bases. As in 0, we use (|0), |1)) to represent basis 1, and 
(|0) = |0) + |1), |1) = |0) — |1)) to represent basis 2. Note that the normalisation 
coefficients have been omitted for the sake of simplicity in presentation. 

A quantum state, denoted by \w), can be coded as 

\ Cw ) = (-l) w » in basis 1, and (3) 

veCi 

\c w ) = \v + w) in basis 2 (4) 

where w G C 2 jC^ \ C 2 jC^ denotes the cosets of in C 2 . Since all the w's 
in Ci define the same quantum state, we can choose w to be coset leaders in 
C^/CiS denoted by id 6 [C^/C^]. There are [C^/C^l = 2 dim ^~ dim ^ coset 
leaders. The quantum code C spans a 2 fc -dimensional Hilbert space where k = 
dim(C 2 1 ) — dim(Cj L ). It is equivalent to encoding a block of k qubits into n qubits 
(mapping from the 2 fc -dimensional Hilbert space into a subspace of a 2 n -dimensional 
Hilbert space). 
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3. Reed-Muller Codes 



Reed-Muller codes are self-dual and hence good candidates for constructing quantum 
error correcting codes ||. In this section, a brief introduction to the structure of 
Reed-Muller Codes is given. A thorough treatment of this subject can be found in 
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For each positive integer m and r (0 < r < m), there exists a Reed-Muller code 
of block length n = 2 m . This code, denoted by RM(r,m), is called the r-th order 
Reed-Muller code. The generator matrix of RM(r, m) is defined as 



G 



G 



(5) 



where Go = {1} is the all-one vector of length n; Gi, an m by 2 m matrix, has 
each binary m-tuple appearing once as a column; and G; is constructed from Gi by 
taking its rows to be all possible products of rows of Gi, I rows of Gi to a product. 
For definiteness, we take the leftmost column of Gi to be all zeros, the rightmost 
to be all ones, and the others to be the binary m-tuples in increasing order, with 
the low-order bit in the bottom row. Because there are C(m, I) ways to choose the 
I rows in a product, G/ is a C(m,l) by 2 m matrix. For an r-th order Reed-Muller 
code, the dimension of the code is given by k = Z)[=o C(m, I). 

Equation fl5|) shows that RM(r— 1, m) is generated by [G , G 1; • • • , G r _i] T , therefore 
RM(r — l,m) C RM(r,m). More generally, 



RM(r -i,m) C RM(r, m) 



for integers 



1 < i < r. 



(6) 



The self-duality property and the minimum Hamming distance of a Reed-Muller 
code can be easily derived from the squaring structure of the code. According to 



Forney jlql , any r-th order Reed-Muller code of length 2 m can be generated through 



recursive squaring construction: 

RM(r, m) = \RM(r, m - 1) /RM{r - 1, m - 1) | 2 (7) 
or two-level squaring construction: 

RM(r, m) = \RM(r, m - 2) /RM(r - 1, m - 2) /RM{r - 2, m - 2) | 4 (8) 
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where RM(r, m - 1) /RM(r - 1, m - 1) and i?M(r, m - 2) /RM(r - 1, m - 2) /RM(r - 
2,m — 2) represent a one-level partition and a two-level partition, respectively. The 
boundary conditions are: 

flMM) = ((M) »r>0 

I (1,0) otherwise 
The squaring construction of RM(r, m) is defined as 

\RM(r,m- l)/RM(r - l,m - 1)| 2 
= {(ti + c,t 2 + c) : 

t h t 2 e RM(r - 1, m - 1), c e [i?M(r, m - l)/RM(r - 1, m - 1)]}. (9) 

From this construction, it is obvious that the minimum Hamming distance of RM(r, m), 
denoted d[RM(r,m)], is given by 

d[RM(r,m)} = mm{d[RM(r - 1, m - 1)], 2d[RM(r, m - 1)]}. (10) 

From the boundary condition, one can easily prove by induction that the minimum 
Hamming distance of RM(r,m) is indeed 2 m_r . 

The dual partition chain of RM(r,m — l)/RM(r — l,m — 1) is i?M ± (r — 1, m — 
1) / i?M _L (r, m — 1). The squaring construction of this dual partition chain is written 

as 

\RM ± (r - l,m- l)/ J RM ± (r,m - 1)| 2 

= {(t' 1 + c ',4 + c ') : 

e RM L {r,m- l),c' e [i?M x (r - l,m - 1) / 'RM ± (r,m - 1)]}. (11) 

| i?M x (r - 1, m - 1) /RM 1 - (r, m - 1) | 2 is the dual of | RM(r, m-1) /RM(r - 1 , m - 1) | 2 
because the inner product of the vectors from the two constructions is zero, 

(*i + c, t 2 + c) • (*i + c', 4 + c') = (c, c) • (c', c') = 0. (12) 

If RM L (r, m-1) = RM[(m - 1) - r - 1, m - 1) for r < m - 1, Eqs (0) and flU) 
shows that RM^ir^m) = RM(m — r — l,m) for r < m. Using the induction method 
and the boundary condition for the initial partition, it can be shown that the dual 
code of RM(r, m) is RM(m — r — 1, m). If m — r — 1 < r, 

RM 1 - = RM(m - r - 1, m) C RM(r, m). 

That is, Reed-Muller codes are self-dual. The minimum Hamming distance of the 
dual code RM(m — r — 1, m) is 2 r+1 . Reed-Muller Codes of block length 4 to 1024 
are listed in Table [L| 
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4. Construction 

Since Reed-Muller codes are self-dual, we have C 1 - C C if k L < k. Let G\ = 
(n,k,d) = (2 m , Z)[ =0 C(m, /), 2 m_r ) and > |. To construct a quantum code from 
Reed-Muller codes, we simply choose C 2 = C±. Therefore, C 2 = (n, fc- 1 , c?- 1 ) = 
(2 m , Y,h=o r ~ 1 C(m, I), 2 r+1 ) and d 1 - > d. Applying the method introduced in Section 
2 and Equations we will have the quantum Reed-Muller code C defined as 

r m—r— I 

C = [[n,k,d}] = [[2 m ,J2C(rn,l)- £ C(m,l),2 m - r ]] (13) 

1=0 1=0 

For example, let m = 10 and n = 2 m = 1024. There are 10 classical Reed-Muller 
codes of length 1024: 

(1024,1023,2) (1024,1013,4) (1024,968,8) (1024,848,16) (1024,638,32) 
(1024,1,1024) (1024,11,512) (1024,56,256) (1024,176,128) (1024,386,64) 

Every pair of dual codes form a quantum Reed-Muller code: 

(1024,1022,2) (1024,1002,4) (1024,912,8) (1024,772,16) (1024,252,32) 

The (1024,252,32) code can correct 15 errors out of 1024 qubits. A comparable 
code in |J is the (1024,462,24) code, which is about 1.8 times more efficient than 
the (1024, 252, 32) code. But the (1024, 252, 32) code is able to correct 4 more errors 
than the (1024,462,24) code. 

A list of new quantum Reed-Muller codes of length 4 to 1024 is given in Table |2|. 
These codes, together with codes listed in ||, form one family of quantum Reed- 
Muller codes. 

Assume that the decoherence process affects each qubit independently, and that 
the error probability of uncoded qubits is p. Then the probability of each coded 
quantum state in error shall be bounded by 

P e < E C(n,j)p n - j (l- P y. (14) 
j=t+l 

The probability of each qubit being in error is given by 

P, = l-(1-P e )» (15) 

The error performances of various quantum error correcting codes are illustrated in 
Figs.|l] and H by applying Equations ( |i~4"D and (|l^). It was found that all single error 
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correcting quantum codes have similar qubit error performance (close to that of the 
[[5, 1, 3]] code). The most effective repetition codes were proposed by Calderbank et. 
al. [§]. Fig.|TJ shows that the [[1024, 252, 32]] quantum Reed-Muller code outperforms 
the [[13, 1, 5]] repetition code significantly, and outperforms the [[29, 1, 11]] repetition 
code asymptotically. The [[1024, 252, 32]] code is about 3 times more efficient than 
the [[13, 1, 5]] code, and 7 times more efficient than the [[13, 1, 5]] code. It should be 
noted though that codes such as [[1024, 252, 32]] are more complex to decode. 

The quantum Reed-Muller codes constructed in this paper are a factor of 2 less 
efficient than those constructed by Steane ||. However they are two orders of 
magnitude more effective, as shown in Fig.|2|. As indicated in both figures, the 
average qubit error probability can be reduced to less than 10~ 9 (one in a billion) if 
the uncoded qubit error rate is not more than 0.3%. 
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Table 1: Parameters of the Classical Reed-Muller Codes 



k 






d 








2 


4 


8 


16 


32 


4 
8 


2 
6 











16 


14 


6 








n 32 


30 


20 









64 


62 


50 


20 






128 


126 


118 


68 







256 


254 


238 


184 


70 




512 


510 


492 


420 


252 





1024 


1022 


1002 


912 


772 


252 



Table 2: Parameters of the Quantum Reed-Muller Codes 
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Figure 1: Qubit Error Performance of repetition codes and the (1024,252,32) Reed- 
Muller code. 
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Figure 2: Qubit Error Performance of Multiple Error Correcting Codes. 
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